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Global aerosol optical models and lookup tables for the new MODIS aerosol 

retrieval over land 

Robert C. Levy, Lorraine A. Remer, Oleg Dubovik 

Since 2000, MODIS has been deriving aerosol properties over land from MODIS 
observed spectral reflectance, by matching the observed reflectance with that simulated 
for selected aerosol optical models, aerosol loadings, wavelengths and geometrical 
conditions (that are contained in a lookup table or ‘LUT’). Validation exercises have 
showed that MODIS tends to under-predict aerosol optical depth (x) in cases of large x (x 
>1.0), signaling errors in the assumed aerosol optical properties. Using the climatology of 
almucantur retrievals from the hundreds of global AERONET suhphotometer sites, we 
found that three spherical-derived models (describing fine-sized dominated aerosol), and 
one spheroid-derived model (describing coarse-sized dominated aerosol, presumably 
dust) generally described the range of observed global aerosol properties. The fine- 
dominated models were separated mainly by their single scattering albedo (co 0 ), ranging 
from non-absorbing aerosol (co 0 ~0.95) in developed urban/industrial regions, to neutrally 
absorbing aerosol (co 0 ~0.90) in forest fire burning and developing industrial regions, to 
absorbing aerosol (co 0 ~0.85) in regions of savanna/grassland burning. We determined the 
dominant model type in each region and season, to create a 1° x 1° grid of assumed 
aerosol type. We used vector radiative transfer code to create a new LUT, simulating the 
four aerosol models, in four MODIS channels. Independent AERONET observations of 
spectral x agree with the new models, indicating that the new models are suitable for use 
by the MODIS aerosol retrieval. 
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1. Introduction 


One of the great scientific uncertainties is the role of aerosols within earth’s 
climate system (IPCC, 2001). Aerosols are deeply involved in the radiation budget, cloud 
processes and air quality. Satellites are increasingly being used for observing global 
aerosol properties. Since MODIS’ launch aboard Terra (in late 1999) and aboard Aqua 
(in early 2002), MODIS spectral reflectance observations have led to retrievals of aerosol 
optical depth (AOD or t) and a measure of the aerosol size distribution, known as the fine 
model weighting (FW or rj). Indeed, the most comprehensive aerosol dataset over land 
has been provided by MODIS (Remer et al., 2005), and such data has been used in 
dozens of applications and publications since launch. Not only have MODIS aerosol 
products been used to answer scientific questions about radiation and climate (e.g. IPCC, 
2001; Yu et al., 2005), they are being used for applications such as monitoring surface air 
quality for health (e.g. Chu et al., 2003, Al-Saadi et al., 2005). 

The recent operational version over land (V4.2) and the products created for 
Collection 004 (C004) were described in Remer et al., (2005). Even though the MODIS 
C004 aerosol products have been ‘validated’ by comparison with sunphotometer, they 
still show room for improvement, especially in certain environments. For example, Levy 
et al., (2005) determined that when compared to sunphotometer data, retrievals of AOD 
over the U.S. East Coast tended to over-predict x in clean conditions, and under-predict in 
more hazy conditions. They postulated that the over-estimation in clean conditions was 
a result of poor land surface reflectance assumptions, and that the under-estimation of 
high x was related to poor assumptions of aerosol optical properties. In a companion 
paper, Levy et al., (2006) discuss surface reflectance assumptions and introduce a new 
retrieval philosophy from MODIS. In this paper, we concentrate on improving the aerosol 
model assumptions. 

The upward spectral ‘reflectance’ (normalized solar radiance) observed by a 
satellite at the top of the atmosphere (TOA) is a function of successive orders of radiation 
interactions within the coupled surface-atmosphere system. The observed spectral 
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reflectance results a combination of processes, including: scattering of radiation within 
the atmosphere without interaction with the surface (known as the ‘atmospheric path 
reflectance’), the reflection of radiation off the surface that is directly transmitted to the 
TOA (the ‘surface function’), and the reflection of radiation from outside the sensor’s 
field of view (the ‘environment function’). The environment function is neglected so that 
to a good approximation, the angular (function of solar zenith, sensor zenith and 
sun/sensor relative azimuth) TOA reflectance at a wavelength A is described by: 


pl(9 0 A<P) = pI(9 0 ,9,<P) + 


s xPx{9 0 ,9,4>) 


where F d A is the ‘normalized downward flux’ for zero surface reflectance, T x represents 
‘upward total transmission’ into the satellite field of view, s x is the ‘atmospheric 
backscattering ratio’, and ff x is the angular ‘surface reflectance’. Except for the surface 
reflectance, each term on the right hand side of Equation 1 is a function of the aerosol 
type (chemical composition, size distribution) and its columnar loading (x). 

The MODIS algorithm follows a lookup table (LUT) strategy, thereby assuming 
that a small set of aerosol types, loadings, and geometry can span the range of global 
aerosol conditions. The LUT contains pre-computed simulations of the non-surface terms 
in equation 1, for the assumed set of aerosol and geometrical conditions. The goal of the 
algorithm is to match the lookup table with the MODIS -observed spectral reflectance ff x , 
to retrieve the associated aerosol properties (including x and ti). The difficulty lies in 
making the most appropriate assumptions about both the surface and atmospheric 
contributions. 

As explained by Remer etal., (2005), the MODIS algorithm over land chooses 
from a set of fine-dominated aerosol models and a single coarse-dominated aerosol 
model. The selection of which fine-dominated aerosol model is fixed based on season and 


location. The coarse-dominated model (dust) is considered fixed, globally. For earlier 
versions (before V4.2) of the retrieval algorithm had a choice of two fine models 
(Kaufman et al., 1997), the ‘urban/industrial’, and the ‘biomass burning/developing 
world’, differing as to their refractive indices, single scattering albedos and phase 
functions. Each of these aerosol models was actually comprised of two or more 
lognormal modes (Kaufman et al., 1997), with their optical properties based on a 
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combination of laboratory studies and sunphotometer data (e.g. Remer et al., 1998). 
Ichoku et al., (2003) found that neither of the two fine models was sufficient to simulate 
the highly absorbing smoke seen over southern Africa. As a result, a third fine model 
(‘highly absorbing’) was added for use in Africa and other absorbing aerosol regions 
(Remer et al., 2005). For C004, the addition of the new ‘highly absorbing’ model 
improved the MODIS/sunphotometer comparison in southern Africa. Although never 
implemented into C004, Levy et al., (2005) showed that revising the urban/industrial 
model could help improve retrievals of high x. It is expected, then, that retrievals in other 
regions (and global) would be improved by re-evaluating the aerosol optical models for 
MODIS. 

Chapter 2 summarizes the new aerosol algorithm. Chapter 3 describes the 
AERONET data used to develop the new models. Chapter 4 discusses cluster analysis of 
the AERONET data. Chapter 5 fixes the aerosol type at given locations as a function of 
season. Chapter 6 introduces the radiative transfer and calculations for the new lookup 
tables. 


2. The C005 aerosol retrieval algorithm 

Levy et al., (2006) introduced a new aerosol algorithm (known as Version 5.2 or ‘V5.2’) 
for deriving aerosol properties from MODIS over land, intended for Collection 5 (C005). 
This new algorithm over land uses gas-absorption-corrected reflectance data in four 
spectral channels (0.47, 0.66, 1.24, and 2.12 pm; MODIS channels 3, 1, 5 and 7), retrieve 
total spectral ‘aerosol optical depth’ (AOD or x) and ‘Fine aerosol Weighting’ (FW or T]), 
reported at 0.55 pm, and ‘surface reflectance’ at 2.12 pm. In contrast to the previous 
family of algorithms (for C004 - Remer et al., 2005), ‘V5.2’ does not need to assume that 
aerosol is transparent in the 2.12 pm (i.e. surface reflectance = observed reflectance at 
2.12 pm). Also, it is not necessary to assume that the surface reflectance values, in the 
two visible channels (0.47 and 0.66 pm), are fixed ratios of that at 2.12 pm. 

Instead, ‘V5.2’ assumes that the relationship between the visible and 2.12 pm is a 
function of scattering angle and a Vegetation Index based on the 1.24 and 2.12 pm 
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observations. This means that the algorithm has some sensitivity to the type of angular 
and surface type effects seen by investigators such as Gatebe et al. (2001) and Remer et 
al. (2001). Using the assumed surface reflectance relationship, the algorithm inverts three 
of the channels (0.47, 0.66 and 2.12 pm) to retrieve the three pieces of information (x, r| 
and surface reflectance). 

Like all previous versions of the over-land algorithm, as well as the MODIS over 
ocean algorithm (Remer et al., 2005; Tanre et al., 1997), the new algorithm is based on 
the lookup table (LUT) strategy. This means that the MODIS observations are compared 
with simulations made by radiative transfer (RT) code, for a small set of aerosol and 
angular conditions that ‘represent’ the possible array of global scenarios. The purpose of 
this paper, then, is to describe the aerosol optical models and RT that leads to the new 
LUT that is used in for deriving C005 products. 

3. The AERONET L2A dataset 

The sunphotometers of the Aerosol Robotic NETwork (AERONET- Hoi ben et al., 
1998) provide a comprehensive data set of aerosol properties. Operating at hundreds of 
sites globally, AERONET has been reporting at some sites since 1993 
(http://climate.gsfc.nasa.gov). ‘Sun’ products are retrievals of spectral x at four or more 
wavelengths (0.44, 0.67, 0.87 and 1.02 pm - some instruments provide more) resulting 
from observations of the spectral extinction of the direct sunbeam. They are provided 
every 15 minutes during the daytime. ‘Sky’ products are hourly estimates of aerosol size 
distribution, phase function and absorption properties, derived from the observed 
scattered spectral radiance in the almucantur. The sky products are derived about Once 
per hour, and assume spherical aerosol particles (Dubovik and King, 2000) and/or 
spheroid particles (Dubovik et al., 2002; Dubovik et al., 2006). Under either particle 
assumption, the fundamental derived parameters include spectral x, spectral complex 
refractive index, the volume distribution as a function of 22 radius size bins (dV/dlnR), 
and fitting error to the radiance measurements. Additional parameters are then calculated 
that include Angstrom exponents, properties of two (fine and coarse mode) lognormal 


6 



aerosol distributions, spectral single scattering albedo (SSA or a> 0 ) and asymmetry 
parameter (g) of the lognormal modes. These data go through rigorous calibration and 
cloud screening processes, resulting in ground-truth estimates of a number of aerosol 
properties. These are known as Level 2 AERONET retrievals (in this paper we designate 
them as ‘L2A’). 

4. Deriving new aerosol optical models 

A number have studies (e.g. Chu et al., 2002, Remer et al., 2005, Levy etal., 2005) have 
demonstrated that MODIS/AERONET regression of x over land results in slope less than 
one; meaning that MODIS tends to under-retrieve x, especially as x increases. Ichoku et 
al., (2003) found that an aerosol model with larger absorption (lower (%) was necessary 
to simulate smoke produced by savanna fires over southern Africa. This highly 
absorbing smoke model was later included for deriving C004. Over the east coast of the 
United States, Levy et al. (2005) showed MODIS could be improved by using the 
urban/industrial aerosol model derived over Goddard Space Flight Center (GSFC) from 
AERONET data (Dubovik et al., 2002). These studies, along with other anecdotal 
evidence, suggested that the set of aerosol models should be re-evaluated. Even though 
the current L2A dataset isn’t truly “global,” it is the most comprehensive ground based 
quality-assured dataset available. 

Omar et al., (2005) performed a “cluster analysis” of AERONET data and found 
that six aerosol models (composed of desert dust, biomass burning, background/rural 
polluted continental, marine, and dirty pollution, respectively) are sufficient for 
representing the entire AERONET dataset. They vary mainly by their oo 0 and size 
distribution. Two models are representative of very clean conditions (marine and 
background/rural). One of the models (dust) is coarse-dominated model, analogous to the 
MODIS coarse dust model, and three are fine models having different co 0 (biomass 
burning, polluted continental, and dirty pollution), that are analogous to the C004 set of 
fine models. . 

Because the MODIS over-land retrieval employs only three channels (and suffers 
from surface and other contaminations), it is not able to select among choices of fine 
aerosol model. Therefore, the algorithm must assign the fine aerosol model apriori of the 
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retrieval. Unfortunately, the Omar et al. study (2005) leaves us one step short of this goal, 
since a unique fine model was not found at every site. Therefore, for use in selecting site 
unique aerosol models, we engaged in a “subjective” cluster analysis that would yield us 
designated models at each site. 

4.1 . Subjective Cluster Analysis of AERONET data 

For our subjective cluster analysis, we downloaded about 136,000 L2A 
almucantur retrievals that were processed as of February 2005. These encompassed both 
spherical and spheroid retrievals. We discriminated the retrievals by the minimum quality 
parameters suggested by the AERONET team, including: t at 0.44 pm greater than 0.4, 
solar zenith angle greater than 45°, 21 symmetric left/right azimuth angles, and radiance 
retrieval error less than 4%. The resulting data set was comprised of 13,496 spherical 
retrievals and 5128 spheroid retrievals. 

In order to differentiate between aerosol types, we separated the AERONET data 
set into ten discrete bins of x. Each bin, then, was used separately to differentiate aerosol 
types. Presumably, this would help to identify expected ‘dynamic’ properties (function of 
AOD) of each aerosol type (e.g. Remer et al., (1998)). For clustering, we employed the 
cluster analysis routines provided with the IDL (Interactive Display Language) software 
version 6.1. Beginning with the entire data set (separately each bin of x), we clustered 
with respect to a number of aerosol optical parameters. In contrast to Omar et al. (2005), 
we desired to pursue not necessarily the most statistically significant clustering, but rather 
to identify distinct models useful for MODIS. With the goal of fine model identification 
in mind, we decided on clustering with respect to only two optical parameters: SSA (<o 0 ) 
at 0.67 pm and the asymmetry parameter (ASYM or g) at 0.44 pm. Presumably co 0 
separates non-absorbing aerosols (such as urban/industrial pollution - (Remer et al., 

1998; Dubovik et al., 2002)) from much more absorbing aerosols (such as savanna 
burning smoke - (e.g. Ichoku et al., 2003; Dubovik et al., 2002)), and g at 0.44 pm would 
help differentiate between the phase functions of different (mainly fine mode - size 
similar to wavelength) aerosols. In the spirit of the C004 aerosol map (Remer et al., 

2005), we looked for three clusters that represented three fine aerosol models. 
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Assuming three clusters for the entire AERONET dataset led to two clusters with 
thousands of points each and one with only 11 points. Presuming the eleven-point cluster 
to contain only outliers, we removed these points from the dataset, and also from the 
appropriate x bins. Because the outliers were removed, now the cluster routine would not 
have to reach so far to find a place for these points, and could be performed again. Re- 
clustering resulted in three clusters for each x bin, each having a reasonable number of 
points. We assumed that we were allowed to number the clusters so that in each x bin, 
one cluster represents the combination of highest co 0 and highest g (a ‘non-absorbing’ 
aerosol model), one cluster represents the lowest co 0 and lowest g (an ‘absorbing’ aerosol 
model), and the third represents the middle combination values (the ‘neutral’ aerosol 
model). As for the coarse aerosol model, we found that a single cluster described the 
spheroid-based almucantur inversions (Dubovik et al., 2006). Since the sites contributing 
to spheroid data were primarily those known to be in dust regions, we assumed that the 
spheroid model represented coarse-dominated (presumably dust) aerosol. 

4.2. Aerosol models around the globe 

Like the C004 algorithm, the C005 algorithm must decide which of the three fine- 
dominated aerosol type to assume for the retrieval. Remer et al., (2005) showed how each 
aerosol type was designated into regions, including places where little was known about 
the prevailing aerosol type. For example, C004 assumed that smoke from tree forest fires 
(both tropical and high latitude forest) had the same optical properties as eastern Europe, 
most of Asia, and other developing regions. How should assumed aerosol type be 
distributed for C005? 

The first step is to determine the aerosol type that best represents each AERONET 
site. For each site, and for each season, we computed the percentage of the retrievals 
attributed to each cluster. Figure 1 (a-d) displays pie-plots at each site, as a function of 
season. To remove poor statistics, we show pie plots only at sites having at least 10 
observations (per season) during the history of AERONET. Unfortunately, this removes 
the many sites that have few retrievals of x > 0.4. Green pie segments represent the non- 
absorbing co 0 ~0.95 model (presumably urban/industrial aerosol), blue segments are the 
neutral oo 0 ~0.90 model (presumably generic, forest smoke and developing world aerosol), 
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and red segments designate the highly absorbing co 0 ~0.85 model (presumably 
savanna/grassland smoke aerosol). At most sites and most seasons, the aerosol type is as 
expected. Non-absorbing aerosol (green) dominates the U.S. East Coast and far western 
Europe, whereas highly absorbing aerosol (red) dominates the savannas of South 
America and Africa. Most other sites are dominated by neutral aerosol (blue) or are a mix 
of all clusters. 

There are some exceptions to expectation, however. Surprisingly, Southeast Asia 
seems to be primarily non-absorbing aerosols, as opposed to the absorbing aerosol 
assumed for C004. Recent studies (e.g. Eck et al., 2005) confirm that aerosol in urban 
areas in far Southeast Asia is primarily non-absorbing. A few sites in Western Europe 
have large fractions of absorbing aerosol, yet the reason is not known. 

Keeping in mind our goal of dividing the world into plausible aerosol types, we 
decided that each site should have an assumed aerosol type attached to it The Neutral 
aerosol model was set as the default, and would be changed only if clear dominance of 
one of the other two aerosol types was observed. If either the non-absorbing or the 
absorbing aerosol occupied more than 40% of the pie, and the other occupied less than 
20%, then the site was designated as the dominant aerosol type. For example, GSFC 
(39°N,77°W) during the summer months (JJA), recorded 87% non-absorbing and 13% 
neutral, meaning it would be designated as non-absorbing. * 

Figure 2 (a-d) displays the designated aerosol types at each site. As in Fig. 1, 
green represents non-absorbing, blue represents neutral and red designates absorbing 
aerosol types. Most site designations seem reasonable and were expected from our 
experience. North America during the summer (JJA) is split between non-absorbing and 
neutral aerosol types, much the same way (approximately -100° longitude) as was 
prescribed for C004. Southern Africa during the winter season (DJF) is solidly 
designated as absorbing aerosol (e.g. Ichoku et al., 2003). Western Europe is evenly split 
between non-absorbing and neutral (except for two absorbing sites), meaning that a 
subjective decision is needed here. To follow the C004 lineage, we chose the non- 
absorbing aerosol model for Western Europe. 

Figure 3 plots the final decision for designating aerosol types around the globe, as 
a function of season. Note that where possible the shapes correspond with the clustering. 
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At some regions, however, some subjectivity was needed to connect areas. For example, 
even though insufficient data exists for Africa north of the equator, the known surface 
types and seasonal cycles suggest that heavy absorbing aerosol would be produced during 
the biomass burning season. Red designates regions where the absorbing aerosol is 
chosen, whereas green represents non-absorbing aerosol. The neutral (co 0 ~ 0.90) model is 
assumed everywhere else. These images were mapped onto a 1° longitude x 1° latitude 
grid, such that a fine aerosol type is assumed for each grid point, globally. This global 
map approach, that is not hardwired into the processing code, will allow for easy 
alterations as new information becomes available. 

Table 1 displays the optical properties and size distributions for the Continental 
model, the three spherical (neutral, absorbing and non-absorbing) fine models and the one 
spheroid (e.g. Dubovik et al., 2006) coarse aerosol (dust) models. Figure 4 shows the size 
distributions for the four AERONET -derived models. Note the dynamic nature (function 
of x) of the size properties of the fine models, especially the non-absorbing model. 

Figure 5 plots the spectral dependence of x (Fig 5A) and phase function at 0.55 
jim (Fig 5B) for each model having x 055 = 0.5. Note that even though the three fine- 
dominated models have similar x spectral dependence, their phase function differs 
somewhat. The coarse model (spheroid-dust) has much smaller spectral dependence than 
any of the fine-dominated models, and nearly flat phase function in the 90°-180° 
scattering angle range observable by MQDIS. 

Figure 6 compares the phase function of each of the new C005 models as 
compared to the analogous models from the C004 algorithm. Changes are minimal 
(especially for the 90°-180° scattering angle range) for the non-absorbing (urban) and 
absorbing (heavy smoke) aerosol types. A possibly significant change is seen in the 
neutral (generic/moderate smoke) phase function. The most obvious change is in the 
‘dust’ models, due to assuming spheroids instead of spheres. The differences are 
primarily in the MODIS scattering angle ranges, which could have a significant effect 
within the aerosol retrieval. 
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5. Radiative transfer and new lookup tables 

5.1 . Channel center wavelengths and Rayleigh optical depths 

The MODIS 0.47 pm blue band (channel 3) stretches between 0.459 and 0.479 
[Am (http://modis.gsfc.nasa.gov). The sea-level Rayleigh optical depth {ROD or x R ) drops 
drastically over this channel, from about 0.203 at 0.459 pm to 0.170 at 0.479 pm 
(Bodhaine et al., 1999). The choice of ‘center’ wavelength to model and its associated 
ROD is crucial to obtaining unbiased aerosol retrieval. This is also an issue for the red 
0.66 pm channel (channel 1: 0.620 - 0.670 pm), but since the RODs are only about one- 
quarter as in the blue, the decision is much less crucial. In C004, the assumed ROD was 
0.186 for channel 3 and 0.048 for channel 1. The 6S RT code (Vermote et al., 1997) 
models the MODIS channel filter functions, and suggests that the ROD values should be 
more like 0.193 and 0.051, respectively for the two channels. The MODIS aerosol over 
ocean algorithm (Taiire et al., 1997, Ahmed et al., 1981) assumes the RODs for the 
channels as 0.195 and 0.052, respectively. 

For C005, re-evaluation of the MODIS channel filter functions showed that 
central wavelengths for channels 3 and 1 are 0.466 and 0.644 pm respectively. According 
to Bodhaine et al. (1999), associated sea-level ROD values are 0.194 and 0.052, 
respectively, which leads to consistency with the aerosol over ocean algorithm. Similar 
calculations for the 0.55 and 2.12 pm channels (channels 4 and 7), found center 
wavelengths of 0.553 and 2.119 pm and RODs of 0.092 and 0.0004, respectively. Note 
that although the center wavelengths are known, we will continue to designate MODIS 
channels 3, 4, 1 and 7 as the 0.47, 0.55, 0.66 and 2.12 pm channels, for brevity and 
consistency with common usage. Table 2 lists the MODIS aerosol channels, along with 
the Rayleigh optical depth assumed for the band. 
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5.2. Radiative transfer code 

The C004 (and previous) LUTs were calculated using ‘SPD’, the scalar version of the 
RT code written by Dave et al., (1970), a code that is a standard in the remote sensing 
community. However, Fraser et al (1989) suggested and Levy et al., (2004) demonstrated 
that under some geometries, neglecting polarization would lead to significant errors in 
top of atmosphere reflectance, further leading to significant errors (> 10% or > 0. 1) in x 
retrieval. Dave also provided a vector (polarized) option to the code (VPD), although it 
has not been well maintained within the MODIS community. Therefore, we need a 
vector code that is well understood and suitable for creating the LUT. In scalar mode, our 
choice of RT code should be consistent with the Dave benchmark. Also, it should 
reasonably match the Ahmed et al., (1981) calculations used for the over ocean aerosol 
retrieval (Remer et al., 2005). 

Levy et al., (2004) employed RT3, the polarized radiative transfer model of Evans 
and Stephens (1991), This plane-parallel adding/doubling code allows for polarization to 
be turned on or off, by changing only one line within an input file. Thus, it was easy to 
compare the results to the Dave code’s scalar mode, and then upgrade to vector mode to 
include polarization effects. Under most geometries and optical depths, differences 
between the two RT codes are less than 0.001 (which is about 1%). 

As in Levy et al., (2004), the spherical aerosol scattering phase function elements 
(inputs to RT3) were calculated by integrating (over size distribution) the results of 
MIEV Mie code (Wiscombe et al., 1980) for hundreds of discrete radii. Properties of the 
aerosol size distribution and refractive indices were those described by Table 1. For the 
spheroids of the coarse aerosol model, Mie theory is not sufficient. We used instead, a 
version of the T-matrix code described in Dubovik et al., (2002, 2006), to calculate the 
scattering properties of the aerosol model. Not only is this a necessary approximation for 
integrating a spheroid size distribution, it is consistent with the calculations used in fitting 
the original almucantur radiance in the first place. 
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5.3. Description of the lookup table (LUT) 

As introduced by Levy et al., (2006), the new V5.2 algorithm over land performs a 
simultaneous inversion of three channels (0.47, 0.66 and 2. 12 pm) to retrieve x, rj and the 
surface reflectance. The inversion technique requires that the LUT be ‘indexed’ like for 
the over-ocean algorithm (Tanre et al., 1997; Remer et al., 2005). Since the principal 
product is x at 0.55 pm, the LUT is indexed in relation to this channel. 

Therefore, the C005 LUT is computed at the four central wavelengths (0.466, 0.553, 
0.644 and 2.119 pm) representing the MODIS channels 3, 4, 1 and 7. The aerosol model- 
dependent parameters of equation 1 are calculated for several values of aerosol total 
loadings (indexed by rat 0.55 pm), and for a variety of geometry. Each of the spherical 
aerosol models (Continental, neutral, absorbing and non-absorbing) and the one spheroid 
model (dust) are represented within the LUT. 

The scattering/extinction efficiencies (Q sca and QJ) of the aerosol size distributions 
are calculated by either MIEV or the Dubovik T-matrix code, depending on the assumed 
shape. Assuming a Rayleigh atmosphere and realistic layering of the aerosol, the 
Legendre moments of the combined Rayleigh/aerosol are computed for each layer of a 
US Standard Atmosphere (U.S. Government, 1976). These moments are fed into RT3 to 
simulate TOA reflectance and total fluxes. 

The scattering and reflectance parameters are calculated for seven aerosol 
loadings (t 055 = 0.0, 0.25, 0.5, 1.0, 2.0 3.0, and 5.0). TOA reflectance is calculated for 9 
solar zenith angles (9 0 = 0.0, 6.0, 12.0, 24.0, 35.2, 48.0, 54.0, 60,0 and 66.0), 16 sensor 
zenith angles (0 = 0.0 to 66.0, increments of 6.0), and 16 relative azimuth angles (<p = 0.0 
to 180.0 increments of 12.0). All of these parameters are calculated assuming a surface 
reflectance of zero. 

When surface reflectance is present, the second term in Equation 1 is nonzero. 

The flux is a function only of the atmosphere, however, the atmospheric backscattering 
term, s, and the transmission term, T, are functions of both the atmosphere and the 
surface. Therefore, RT3 is run two additional times with distinct positive values of 
surface reflectance. 
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s = (l/ri)(l-(F d Tpi/ip*~p a ))) 

and (2) 

S = (l/p° 2 )(l-(F d Tp s 2 /(p-p a ))) 

Here, we chose values of 0.1 and 0.25 for our surface reflectance. /O^and ff ? These two 
equations can be solved for the two unknowns, s and T. The values of F d , s, and T are 
saved into the LUT, for each x index, wavelength and aerosol model. The scattering and 
extinction coefficients, Q, are saved into the LUT. In addition, a parameter known as the 
Mass Concentration coefficient (1VQ is reported (see Appendix for its derivation). 

6. Spectral dependence of aerosol models 

How well do our derived models represent ambient aerosol at specific AERONET 
sites? For this purpose, we used the time series of Level 2 ‘sun’ retrieved products from 
AERONET, which are independent of the ‘sky’ retrieved products. We cannot use the 
sun measurements to evaluate the assumed absorption properties, but we can analyze the 
spectral dependence of the aerosol optical depth. 

The LUT is indexed by optical depth at 0.55 pm, in increments including x = 

0.25, 0.5, 1.0 and 2.0. Figure 5A showed the spectral x dependence of each model for 
x 0 .55=0.5. Similar plots could be made from the spectral dependence of the other indices 
within the LUT. For each AERONET site, we separated retrievals into three-month 
seasons (winter = DJF, spring = MAM, summer = JJA, fall = SON). Observations within 
each season were sorted according to x 055 , where x 0 55 was calculated by fitting a quadratic 
to the observed spectral x. For each indexed value of x 055 (x = 0.25, 0.5, 1.0, 2.0), we 
determined which AERONET observation contained x^ closest in magnitude to the 
indexed value. This we considered this location the ‘central’ (Q) index of the bin. The 
total number of observations ‘A,’ in each season, divided by 20, determined the number 
i N' of observations that should be considered close to the ‘central’ index. Therefore, the 
set of AERONET observations for the x 0 55 bin spanned between ‘Q - N/2’ and ‘Q + N/2’. 
If there were not enough observations near a certain x 0 55 value, we tried N=N/40. The 
spectral optical thickness for each bin was calculated by averaging the spectral optical 
thickness for the set of observations within the bin. 
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Figure 7 compares spectral dependence of the aerosol models with spectral 
dependence at selected AERONET sites, for selected indexed values of t 055 . The MODIS 
LOT is calculated at four wavelengths (0.47, 0.55, 0.66 and 2.12 pm) channels, whereas 
AERONET spans four to eight bands between 0.34 and 1.02 pm (depending on which 
site). Though not contained in the lookup table, x at 0.855 pm is also modeled for 
MODIS. Interpolation of AERONET to 2.12 pm was not performed because of the great, 
distance from 1.02 pm. Different seasons (for AERONET) are represented by different 
line styles, whereas different t 0 55 indices (values) are color-coded. Because at some sites, 
the assumed aerosol type changes by season, spectral dependence for an additional model 
is plotted. 

At AltaJFloresta (9°S, 56° W), the spectral dependence in the visible wavelengths 
agrees with either the neutral or absorbing models. At 0.855 pm, the AERONET spectral 
dependence varies as a function of season, which somewhat mirrors the differences in 
modeled spectral dependence. During the summer and fall, the AERONET dependence is 
slightly closer to the absorbing model than the neutral, and during the winter and spring, 
the neutral model provides a slightly better match. 

At Capo_Verde (16°N, 22°W), although the neutral fine-dominated model is 
assumed all year, coarse (dust) is expected to dominate. Plotted for Capo_Verde is the 
AERONET spectral t compared with the modeled dust. Even though we believe that we 
have improved the dust model from C004, the modeled spectral dependence is still too 
large to properly represent dust over Capo_Verde. 

The non-absorbing model (co 0 ~ 0.95) shows remarkable match to observations at 
GSFC (39°N, 77° W). The only difference is seen during the winter and spring for the 
lowest x value (0.25), where the particles are known to be larger (have less spectral 
dependence) than the rest of the year. Mongu (15°S, 23°E) is another site that is well 
represented by its assumed aerosol type (absorbing). 

Beijing (39°N, 1 16°E) and Venise (45°N, 12°E) are interesting because dominant 
aerosol type is known to vary. Both sites are influenced by dust transport, so that the 
averaged AERONET spectral dependence should lie somewhere between the fine- 
dominated and coarse-dominated (dust) models. Due to plotting constrains, the coarse 
model (dust) spectral dependence is not plotted, but it is obvious that over Beijing the 
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aerosol is mixed neutral and dust, and more coarse-dominated during the winter and 
spring. Venise is less often in the path of dust (from Africa) but its averaged spectral 
dependence shows the addition of coarse aerosol not represented by a fine-dominated 
model. 

7. Conclusion 

Since 2000, MODIS has been capably retrieving aerosol properties over land. 
However, the long-term study of MODIS products led us to re-evaluate the assumed 
aerosol model physical and optical properties. We used the entire time series of 
almucantur-derived aerosol properties from all AERONET sites to compose a set of three 
fine-dominated (spherical) and one coarse-dominated (spheroid) aerosol optical models 
that seemed to best represent the global aerosol system. The fine-dominated aerosol types 
differed mainly by their values of oo 0 , so se designated them as ‘neutral' (co 0 ~0.90), 
‘absorbing’ (co 0 ~0.85) and ‘non-absorbing’ (oo 0 ~0.95). We then created seasonal l°xl° 
maps of designated aerosol type. 

Phase functions for some of the models were nearly identical to the analogous 
models from C004 algorithm. One of the models (coarse) had substantially different 
phase function from the C004 due to the assumption of spheroids instead of spheres. 
Spectral dependence of x of the models was compared with ‘sun’ retrieved observations 
from selected AERONET sites. Sites dominated by fine aerosol (e.g. GSFC and Mongu) 
were well represented by their assumed aerosol type. At Alta_Floresta, also dominated by 
fine aerosol, our seasonal choice of fine models (neutral versus absorbing) was correct In 
a dust-dominated site such as Capo_Verde, our dust model has too much spectral 
dependence. Sites that are influenced' by occasional dust episodes (e.g. Venise or Beijing) 
show spectral dependence lying between the fine and coarse -dominated models. 

While the derivation of the new set of aerosol models is an interesting exercise in 
its own right, the real strength to the work will be its implementation within the new 
aerosol algorithm (Levy et al., 2006). Preliminary tests of the combination of updated 
aerosol models and new algorithm procedure is expected to improve the retrieved 
products for C005 (Levy et al., 2006). 
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However, we note that the MODIS evaluation is never complete. The AERONET 
data team has nearly completed revision of the sky products dataset. Reprocessing of this 
dataset may lead to updated aerosol models, especially for dust-dominated sites. We also 
plan to evaluate assumptions of oj 0 > using independent datasets. 
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Appendix 1 : Derivation of Extinction, Scattering and 
Mass Coefficients 


The following equations lead to derivation of Mass Concentration in units of (pg per 
cm 2 ). In these equations: dNldXnr is the number size distribution with r denoting radius 
(in pm). For a lognormal mode, r g is the geometric mean radius, o is In a g representing 
the standard deviation of the radius, and N 0 is the number of particles per cross section of 
the atmospheric column (i.e. the amplitude of the lognormal size distribution). In our 
case, we assume that the distribution is normalized, so that N 0 ~\, 


The number N is related to the volume V and area A distributions by: 

dN. 3 _ 3 dV _2 dA 

= —nr — — = Jtr 


dlnr 4 dlnr dlnr 


N 0 , V 0 , and A 0 are the amplitudes of the corresponding distributions, i.e. 

, dV j r* dN dA 

Vo = I dlnr N n = I -dlnr A n = I <ilnr. 

0 Jo dlnr 0 Jo dlnr 0 Jo dlnr 

For a single lognormal mode defined by: 


dN 1 N n 


exp 


/ ln(r/r) 2X 


dlnr r o-yfliz 
and 
3 3 

N 0 = V 0 -Jtr g ex 


2a 2 


) . 




the Moments of order k, Vf are defined as 

M k = f r k — din r = (r ) k exp(0.5 k 2 o 2 ) . 
Jo dlnr 8 


The effective radius r eff in (pm] is defined by the moments, i.e. 
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r eff = 


M_ 

M 2 


f) 

Jo 

T 

Jo 


dN 

dlnr 


dlnr 


* 2 dN „ 
r — dlnr 


3 Vo .5 2 , 

-= r exp(— CT ) 

4 A 0 * F 2 


dlnr 


The extinction coefficient, is related to the extinction efficiency <2 m through the area 
distribution, and is specific to each mode 




r 

Jo 


2 diV 

?ir dlnr 

dlnr 


These parameters are calculated via Mie code (MIEV, (Wiscombe, 1980)). Note that the 
scattering coefficient /3 sca and efficiency Q sca are related the same way. The mass 
extinction coefficient B^, is in units of (area per mass) and depends on the extinction 
efficiency and the particle density p (assumed to be 1 g per cm 3 ), such that 


B ext = 


4pr 


eff 


(Chin et al., 2002). For a single lognormal mode, 

^Qext 3 Qex, A 0 Qext A 0 


B pr , = 

ext 


& 


30. 


4pr e/f 4 p (3/4)V 0 p F 0 pF 0 4 par exp(4.5<7 ) 


However our aerosol models are sums of multiple modes, so that the area and 
volume distributions must take into account the contributions of each mode. If there are 
two modes, (i.e. modes 1 and 2), r e// must be calculated this way: 

r r >m±m d]nr 

Jo dlnr 




r r ,m±M£ dlBr 

Jo dlnr 


Similar modifications are made when calculating Q and thus B. 

We can define the Mass Concentration Coefficient M c , as the inverse of B, such 

that 

1 


M =■ 


B 


The columnar mass concentration, M (mass per area) can then be defined as 
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M - tM c = ■— 
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Listed for each model are the individual lognormal modes, and the final SS A at different wavelengths. Listed for each mode are the 
mean radius r v , standard deviation a of the volume distribution, and total volume of the mode, V 0 . The complex refractive index is 
assumed for all wavelengths (0.47, 0.55. 0.66 and 2.1 pm), unless otherwise noted. The Absorbing and Neutral model parameters (r v , 
a and k) are defined for x < 2.0; for x > 2.0, we assume t = 2.0. Likewise, the Non-absorbing and Spheroid model parameters are 
defined for x< 1.6. V 0 (for all models) is defined for all x. 

TABLE 2: CHARACTERISTICS OF MODIS CHANNELS USED IN THE AEROSOL RETRIEVAL 


Band# 

Bandwidth (pm) 

Weighted Central 
Wavelength (pm) 

Resolution (m) 

Rayleigh optical depth 

1 

0 . 620 - 0.670 

0.646 

250 

0.0520 

2 

0 . 841 - 0.876 

0.855 

250 

0.0165 

3 

0 . 459 - 0.479 

0.466 

500 

0.1948 

4 

0.545 - 0.565 

0.553 

500 

0.0963 

5 

1 . 230 - 1.250 

1.243 

500 

0.0037 

6 

1 . 628 - 1.652 

1.632 

500 

0.0012 

7 

2 . 105 - 2.155 

2.119 

500 

0.0004 


Notes: Band #26 (1.38 pm channel) is used for cirrus correction; NeAp corresponds to the sun at zenith (0 - 0°) 
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Figure 2: Final spherical aerosol model type designated at each AERONET site per 
season. Colors represent absorbing (o) 0 ~0.85), neutral (to 0 ~0.90) and non-absorbing 
(o) 0 ~0.95), respectively. 
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Figure 3: Final spherical aerosol model type designated at 1° x 1° gridbox per 
season. Red and green represent absorbing (to 0 ~0.85) or non-absorbing (co # ~0.95) 
models, respectively. Neutral (to 0 ~0.90) is assumed everywhere else. 
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Figure 4: Aerosol size distribution as a function of x (AOD) bin for the three 
spherical (neutral, absorbing and non-absorbing) and spheroid (dust) models 
identified by clustering of AERONET. 
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Figure 5: Spectral x dependence (A) and phase function at 0.55 pm (B) for the 5 
aerosol models in of the V5.2 LUT. For both plots, to ss = 0.5. 
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Figure 6: Comparison of phase function (at 0.55 pm) between V5 (solid curves) and 
analogous V4 (dotted curves) aerosol models. Models are neutral (A), nonabsorbing 
(B), absorbing (C) and dust (D). For all plots, t 055 - 0.5. 
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Figure 7: Comparison of spectral x between V5 models (filled shapes) and averages 
of AERONET ‘sun’ measurements (dotted curves) at selected sites and seasons. The 
colors represent discrete (Indexed=I) values of t 0 55 , such that green=x o 55 =0. 25=1=1, 
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blue= x 0 5 5 =0.5=I=2, red=nr 0 5 5 =1.0=I=3 and black=T o 55 =2.0=I=4. Different curve 
types represent AERONET data during a season (dot = DJF (W), short dash = 
MAM (P), medium dash = JJA (S) and long dash = SON (F». In the legend, ‘N’ is 
the number of AERONET observations contained within each curve. Note that for 
each site, spectral x of more than one aerosol model and/or fewer than all seasons 
and all x indices may be plotted. Sites plotted are Alta_Floresta (A), Capo_Verde 
(B), GSFC (C), Mongu (D), Beijing (E) and Venise (F). 



